SCIENTIFIC 

REPORTS 




OPEN 



SUBJECT AREAS: 
POTASSIUM CHANNELS 
SENSORS AND PROBES 

Received 
28 January 2014 

Accepted 
24 April 2014 

Published 
12 May 2014 



Label-free cell phenotypic profiling 
decodes the composition and signaling 
of an endogenous ATP-sensitive 
potassium channel 



Haiyan Sun'*, Ying Wei', Huayun Deng', QiaojieXiong^f, Min Li^, Joydeep Lahiri' &Ye Fang' 

'Biochemical Technologies, Science and Technology Division, Corning Incorporated, Corning, NY 1 483 1 , United States of 
America, ^The Solomon H. Snyder Department of Neuroscience and High Throughput Biology Center, Johns Hopkins University 
School of Medicine, Baltimore, Maryland 21 205, United States of America. 



Correspondence and 
requests for materials 
should be addressed to 
Y.F. (fangy2@corning. 



* Current address: 
Biodesign Institute, 
Arizona State 
University, AZ 85287, 
USA 

t Current address: 
Cold Spring Harbor 
Laboratory, Cold 
Spring Harbor, NY 
1 1 724, United States 
of America 



Current technologies for studying ion channels are fundamentally limited because of their inabiUty to 
functionally link ion channel activity to cellular pathways. Herein, we report the use of label-free cell 
phenotypic profiling to decode the composition and signaling of an endogenous ATP-sensitive potassium 
ion channel (Katp) in HepG2C3A, a hepatocellular carcinoma cell Une. Label-free cell phenotypic agonist 
profiling showed that pinacidU triggered characteristically similar dynamic mass redistribution (DMR) 
signals in A431, A549, HT29 and HepG2C3A, but not in HepG2 cells. Reverse transcriptase PGR, RNAi 
knockdown, and K^tp blocker profiling showed that the pinacidU DMR is due to the activation of SUR2/ 
Kir6.2 Katp channels in HepG2G3 A cells. Kinase inhibition and RNAi knockdown showed that the pinacidil 
activated K^tp channels trigger signaling through Rho kinase and Janus kinase-3, and cause actin 
remodeling. The results are the first demonstration of a label-free methodology to characterize the 
composition and signaling of an endogenous ATP-sensitive potassium ion channel. 

I Ion channels are pore forming membrane proteins that coordinate diverse and vital functions including 
contraction, neurotransmission, secretion, and cell volume changes. Despite being implicated in a broad 
spectrum of diseases ranging from neurodegenerative disorders to diabetes, only 7% of all marketed drugs 
are targeted against ion channels'"'. Ion channels remain an underexploited target class in large part due to the 
lack of appropriate screening methods. Current technologies such as automated patch clamping, radioactive flux 
and ion-specific fluorescence dye assays are viable options but suffer from a common drawback - they only offer 
information about ion-channel opening/closing and potential upstream effectors'" \ As a result, our understand- 
ing of cellular pathways linked to ion channels remains poor relative to other membrane proteins such as G 
protein-coupled receptors (GPCRs) and receptor tyrosine kinases. This is exemplified by adenosine triphosphate 
(ATP) sensitive potassium (Katp) channels. Katp channels serve as molecular sensors linking membrane excit- 
ability to metabolism'' The Katp channels are activated by interaction with intracellular Mg^* ADP and inhib- 
ited by high level of ATP, thus are sensitive to the energy state of cells'''". Katp channels express in various tissues 
with different molecular compositions including cardiac myocytes, pancreatic |3 cells, smooth muscle cells, and 
neurons''". Cell plasma membrane Katp channel protein complex consists of four small pore-forming inward 
rectifying potassium channel subunits (Kir6.1 or Kir6.2) and four sulfonylurea receptors (SURl, SUR2A, or 
SUR2B) as the regulatory (3 subunits'"'". Besides expression at the cell plasma membrane, the Katp channels are 
also found in the inner membrane of mitochondria'* and nuclear envelope'^. Recent studies support that cardiac 
mitochondria Katp channels contain renal outer medullary potassium channel (ROMK), beside SUR2; however, 
their exact composition and organization remain elusive"""*. Advances in molecular biology and ion channel 
techniques have deepened our understanding about the assembly, expression, gating, structure, and regulation of 
Katp channels. However, little is known about cellular signaling mediated through Katp channels'^ 

Here we report the presence, composition and signaling of a functional Katp channel in HepG2C3 A cells. This 
discovery was made possible by label-free cell phenotypic profiling with the aid of molecular and chemical biology 
tools. Label-free resonant waveguide grating (RWG) biosensor in microplate allows for non-invasive and real- 
time detection of cell phenotypic responses, termed dynamic mass redistribution (DMR), arising from the ligand- 
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Figure 1 | DMR heat map of ion channel Ugands active in five cell lines, A431, A549, HT29, HepG2, and C3A. This heat map was obtained using 
clustering analysis of the DMR profiles of the ion channel ligands in the five cell lines. For each ligand, the real responses at six discrete time points post- 
stimulation (3, 5, 9, 15, 30, 45 min) in each cell line were used for the cluster analysis. Only ligands that led to a DMR greater than 50 pm, or less than 
— 50 pm in at least one cell line were included in this analysis. 



induced receptor activation in native cells^"'^'. The DMR signal 
obtained offers a holistic view of receptor signaling, so it is possible 
to deconvolute the systems cell biology"""^'' and signaling waves^^ of 
different classes of receptors. Taking advantages of the label-free 
biosensor to be highly sensitive and non-invasive, we used it to 
identify the DMR signature of pinacidil, a K^tp channel opener^*'^'^, 
across multiple cell lines, and deconvolute its origin and signaling 
pathways in HepG2C3A cells. We demonstrated that HepG2C3A 
cells express a functional K^tp channel, although its location remains 
unknown. 

Results 

Label-free cell phenotypic profiling of ion channel ligands. To 

identify functional receptors as well as appropriate cell line(s) for 
studying endogenous Katp channels, we first profiled a 
commercially available library consisting of seventy-two ion 
channel ligands, each at 10 |iM, in the five distinct cell lines using 
DMR agonist assay in microplate^'. The negative controls (that is, the 
cell response to the buffer solution containing 0.1% DMSO, equal 
amount to those for all ligands) were also included. Using the 3ct of 
the negative control response, we limited our analysis to active 
ligands that led to a DMR greater than 50 pm or smaller than 
— 50 pm in at least one cell line. Results showed that out of 
seventy-two ligands twenty-two triggered robust DMR in at least 
one cell line (Fig. 1). Similarity analysis with the Ward hierarchical 
clustering algorithm and Euclidean distance metrics^*" showed that 
these active ligands are divergent in their label-free cell phenotypic 
agonistic activity in these cell lines. 

Detailed analysis revealed several interesting findings. First, thap- 
sigargin, A23187, and cyclopiazonic acid shared similarity in their 
cell phenotypic agonistic activity. AU three ligands triggered robust 
positive DMR in A43 1 , A549, and HT29. However, only thapsigargin 
and A23187 were active in HepG2 and C3A cells. The non-selective 
agonistic activity of these ligands across the five cell lines, together 
with their similarity in DMR characteristics to those arising from the 
activation of Gq-coupled receptors in respective cell lines™, suggests 



that their DMR is due to the increasing cytosolic calcium concentra- 
tions induced by these ligands. Thapsigargin is a non-competitive 
inhibitor of sarco/endoplasmic reticulum Ca^* ATPase (SERCA)"", 
while cyclopiazonic acid is a specific inhibitor of SERCA^^. Both 
inhibitors are known to raise cytosolic calcium concentration by 
blocking the ability of the cell to pump calcium into the sarcoplasmic 
and endoplasmic reticula, which leads to depletion of these stores 
and activates plasma membrane calcium channels, allowing an influx 
of calcium into the cytosoP'. A23187 is a Ca^"^ ionophore used exten- 
sively to mimic the effect of many physiological cell stimuli related to 
calcium ion". The less robust response induced by cyclopiazonic acid 
is consistent with its relatively low potency, comparing to thapsigar- 
gin'"'. Due to its low potency, cyclopiazonic acid at 10 [iM may not 
trigger the maximal response. 

Second, NPPB (5-nitro-2-(3-phenylpropylamino)benzoic acid), 
niflumic acid, IAA-94 (R(-f )-methylindazone) and flufenamic acid 
shared similar cell phenotypic pharmacology; aU were active in HT29 
and C3A cells, and to less extent in HepG2 cells. Recently, we had 
showed that the positive DMR of niflumic acid in HT29 is due to the 
activation of GPR35''''. Consistent with its known agonist activity at 
the GPR35'=, NPPB triggered a dose-dependent DMR in HT29 
(Fig. 2a), yielding an apparent logECso of —5.95 ± 0.07 (n = 3) 
(Fig. 2b). It also desensitized the response to 1 |iM zaprinast, a 
known GPR35 agonist, with a loglCso of -5.79 ± 0.05 (n = 3) 
(Fig. 2b). The known GPR35 antagonist CID2745687 dose-depen- 
dently and completely blocked the DMR of 4 [iM NPPB (n = 3) 
(Fig. 2c), with an apparent loglCso of -5.51 ± 0.06 (n = 3) (Fig. 2d). 
Although flufenamic acid was reported to be inactive in human 
GPR35-P-arrestin-2 interaction assays'*", we found that flufenamic 
acid is a DMR biased partial agonist at the GPR35 (Fig. 2e-h). 
Flufenamic acid triggered a dose-dependent DMR in HT29 
(Fig. 2e), yielding an apparent logECso of —5.18 ± 0.03 (n = 3) 
(Fig. 2f). It also desensitized the cells responding to the second stimu- 
lation with 1 |iM zaprinast with a loglCso of -4.95 ± 0.06 (n = 3) 
(Fig. 2f). CID274568 dose-dependently and partially blocked the 
DMR of 32 |iM flufenamic acid with a loglCso of ^5.28 ± 0.05 
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Figure 2 | Agonist activity of two ion channel ligands, NPPB and flufenamic acid at the GPR35 in HT29 cells, (a) The real time DMR dose responses of 
NPPB. (b) The real DMR amplitudes of NPPB as a function of its dose, in comparison with the DMR of 1 |tM zaprinast as a function of NPPB doses, 
(c) The real time DMR of 4 |tM NPPB in the presence of CID2745687 at different doses, (d) The DMR ampUtudes of 4 |tM NPPB as a function of 
CID2745687 dose, (e) The real time DMR dose responses of flufenamic acid, (f) The DMR amplitudes of flufenamic acid as a function of its dose, in 
comparison with the DMR of 1 )tM zaprinast as a function of flufenamic acid dose, (g) The real time DMR of 32 )tM flufenamic acid in the presence of 
CID2745687 at different doses, (h) The DMR amplitudes of 32 |tM flufenamic acid as a function of CID2745687 dose, (i) The (J-arrestin Tango signal as a 
function of compound dose. For the desensitization, the cells were pretreated with respective ligands for 1 hr before stimulation with zaprinast at the fixed 
dose ( 1 |tM). For the antagonist blockage, the cells were pretreated with CID2745687 at different doses for 1 hr before stimulation with respective agonist 
at a fixed dose. Data represents mean ± s.d. (n = 3). 



(Fig. 2g & h). Using GPR35 Tango P-arrestin gene reporter assay, we 
found that NPPB triggered a dose response with a logECso of —4.55 
± 0.04 and a maximal response that is 45 ± 3% (n = 3) of the full 
agonist zaprinast, but flufenamic acid was inactive (Fig. 2i). These 
results suggest that similar to several other GPR35 agonists including 
benserazide''^, tolcapone"* and rosmarinic acid"", flufenamic acid not 
only activates the endogenous GPR35 in HT29, but also activates an 
additional unknown receptor. Of note, we did not determine the 
mechanism accounted for the DMR of IAA-94 in HT29, given that 
the IAA-94 DMR is much smaller than those induced by GPR35 
agonists. 

Third, punozide, SDZ-201106 and ZM226600 all triggered char- 
acteristically similar DMR specifically in A431. ZM226600 is a Kir6 
Katp channel opener with an EC50 of 0.5 iiM*". Pimozide is a dopa- 
mine D2 antagonist and has high affinity to 5-HT7 receptor'". SDZ- 
201106 is a sodium channel opener*^. Given that no common mode 
of action is known among these ligands, it is possible that they activ- 



ate an unknown receptor in A431 cells; this possibility is currently 
under investigation. 

Fourth, pinacidU was found to trigger a robust negative DMR in 
A431, A549, HT29 and C3A, but not in HepG2 cells. However, 
several other K^^tp openers including diazoxide, minoxidil, minox- 
idil sulphate, and PCO-400 in the library, each at 10 |tM, were inact- 
ive in the DMR agonist profiling across the five cell lines. Given that 
Katp openers are classified based on their ability to open the channel 
per se, but the DMR of a ligand reflects the functional consequence of 
receptor-ligand interaction in cells, we postulated that similar to 
GPCR ligands*^, Katp openers may also display functional selectivity 
and the DMR of pinacidil is related to the signaling of endogenous 
Katp channels. Of note, it was reported that C3A expresses Katp 
channels that play a role in cell proliferation'*'', and HepG2 has little 
Katp channel current on the cell surface but the Katp channel 
expression can be induced by transfection of insulin and glucose 
transporter GLUT2'*l 
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Figure 3 | Dose responses of Katp openers in C3A and A549 cells, (a) The 

real-time DMR dose responses of pinacidil in C3A cells, (b) The real DMR 
amplitudes at 50 min post stimulation as a function of ligand dose in C3A 
cells, (c) The real-time DMR dose responses of pinacidil in A549 cells, 
(b) The real DMR amplitudes at 50 min post stimulation as a function of 
ligand dose in A549 cells. Data represents mean ± s.d. (n = 6 for pinacidil, 
and n = 3 for other openers). 

Together, these analyses suggest that the label-free cell phenotypic 
profiling approach is powerful for determining the on-target and off- 
target activity of ion channel ligands. Given the unique signature of 
pinacidil, we herein were primarily focused on the deconvolution of 
the target and pathways activated by pinacidil. 

The pinacidil DMR in C3A cells is due to the activation of 
endogenous Kj^tp channels. To determine the composition of 
endogenous Katp channels, we first characterized the dose 
responses of known Katp openers using DMR agonist assay. 
Results showed that pinacidil triggered a dose response in C3A 
cells (Fig. 3a), yielding a logECso of -4.77 ± 0.05 (n = 6) 
(Fig. 3b). In comparison, both PCO-400 and diazoxide triggered 
much smaller negative DMR with lower potency; their logECsoS 
were -3.99 ± 0.07 and -3.47 ± 0.09 (n = 3), respectively 
(Fig. 3b). In contrast, minoxidil up to 64 [iM was inactive. 
Similarly, pinacidil also triggered a dose response in A549 cells 
(Fig. 3c), yielding a logECso of -4.72 ± 0.05 (n = 6) (Fig. 3d). In 
comparison, PCO-400 triggered a much smaller negative DMR with 
a logECso of —4.16 ± 0.03 (n = 3), while both diazoxide and 
monoxidU were inactive (Fig. 3d; data not shown). The relatively 
low potency of active Katp openers obtained may be related to the 
whole cell measurements, wherein the intracellular ATP/ADP ratio 
has direct effects. Another possibility is that the Katp channels are 
located inside cells, and the effective intracellular concentration of 
these openers is a function of cell uptake and efflux and may be lower 
than those added extracellularly. Of note, in order to better visualize 
the DMR characteristics we used real DMR responses, instead of 
their absolute values, for aU dose response analysis. Nonetheless, 
these results suggest that these openers display distinct potency 
and efficacy to trigger DMR in C3A and A549 cells. 

Second, we performed reverse transcriptase PCR to determine the 
mRNA expression of Katp channels. Results showed that C3A endo- 
genously expresses mRNAs for SUR2B, Kir6.1, and Kir6.2, to less 
extent SUR2A, but not SURl (Fig. 4a). A549 expresses mRNAs for 



SURl, Kir6.1 and Kir6.2, but not SUR2A and SUR2B (Fig. 4b). A431 
expresses mRNAs for SUR2A, Kir6.1 and Kir6.2, but not SURl and 
SUR2B (Fig. 4c). HT29 expresses mRNAs for primarily SUR2B, 
Kir6.1 and Kir6.2, and to less extent SURl and SUR2A (Fig. 4d). 
These results suggest that different cell lines have different expression 
patterns of Katp channels. 

Third, we applied RNAi knockdown to determine the functional 
Katp channels in C3A cells. Results showed that the mock or 
scrambled RNAi transfection had little impact on the pinacidil 
DMR (Fig. 5a). The treatment with three RNAi against SURl had 
little impact on the pinacidil DMR (Fig. 5b), consistent with the 
absence of SURl mRNA in C3A cells. However, knockdown of 
SUR2 with three RNAi all markedly suppressed the pinacidil DMR 
(Fig. 5c). Similar to the mock or scrambled RNAi transfection 
(Fig. 5d), three RNAi knockdown of Kir6.1 had little impact on the 
pinacidil DMR (Fig. 5e). In contrast, the knockdown of Kir6.2 with 
three RNAi all markedly suppressed the pinacidil DMR (Fig. 5f). 
One-way ANOVA analysis suggests that RNAi for SUR2 and 
Kir6.2, but neither SURl or Kir6.1, significantly altered the pinacidil 
DMR (Fig. 5g and h). Given the low expression level of these proteins 
and the moderate efficiency of typical RNAi knockdown in our 
laboratory^"" (also see below), we didnot attempt to use Western blot 
to quantify the knockdown efficiency of respective RNAi. 

Fourth, the effects of different Katp inhibitors on the pinacidil 
DMR were examined and used to further strengthen RNAi knock- 
down results. Results showed that three sulfonylurea blockers includ- 
ing tolazamide, glipizide and tolbutamine all dose-dependently 
blocked the pinacidil DMR, leading to a single loglCso of —3.85 ± 
0.09, 3.94 ± 0.10, and 3.91 ± 0.10 (n = 3), respectively. In contrast, 
glibenclamide up to 250 |.iM only partially attenuated the pinacidil 
DMR, and U-37883A had little impact on the pinacidil DMR (Fig. 5i 
& j). U-37883A is a non-sulfonylurea blocker that has been reported 
to selectively inhibit Kir6.1 containing Katp channels'"", while glib- 
enclamide is a relatively selective SURl inhibitor*^''"'. Furthermore, 
the co-existence of 40 [iM tolazamide caused the right shift of pina- 
cidil's potency (Fig. 5k). Together, these results suggest that the 
pinacidil DMR in C3A is due to the activation of SUR2/Kir6.2 
Katp channels. Of note, we attempted, but failed to conclusively 
determine the exact location of Katp channels in C3A cells, due to 
very low levels of Katp channel proteins. 

Collectively, these pharmacological profiling, RNAi knockdown, 
and RT-PCR results suggest that C3A endogenously expresses a 
functional SUR2/Kir6.2 Katp channel, whose activation by pinacidil 
results in a robust DMR signal. 

Kinase inhibition profiling. The regulation of Katp channels has 
been extensively investigated; in particular, its upstream regulation 
mediated by protein kinase A (PKA) has been well documented"""^^. 
The phosphorylation of ion channels by protein kinases is an 
important mechanism by which membrane excitability is regulated 
by signaling pathways. However, a clear understanding about the 
signaling of Katp is still lacking. Thus, we investigated the impact 
of a library of known kinase inhibitors on the pinacidil DMR in C3A 
cells. Results showed that the majority of kinase inhibitors had little 
impact on the pinacidil DMR (Fig. 6a). However, a subset of 
inhibitors markedly suppressed the DMR of pinacidil. Among 
them, H-89, HA- 1077, H-7, and H-8 are known PKA inhibitors, 
while AG490 is a putative Janus activated kinase (JAK) inhibitor, 
Y27632 is a ROCK inhibitor, and the alkaloid staurosporine is a 
broad spectrum, high affinity kinase inhibitor with highest affinity 
for protein kinase C (0.7 nM). DMR inhibition assay showed that H- 
89 and H-7 dose-dependently and completely blocked the DMR of 
32 |iM pinacidil, leading to a loglCso of -5.32 ± 0.09 and -5.17 ± 
0.10, respectively (Fig. 6b). The sensitivity of the pinacidil DMR to 
these PKA inhibitors were consistent with the previous findings, 
which suggest that PKA is anchored in proximity to Katp channels 



SCIENTIFIC REPORTS | 4 : 4934 | DDI: 1 0. 1 038/srep04934 



4 



bp 

1000 

500 

300 
200 
100 



HepG2C3A 




A549 




A431 



bp 

1000 

500 
300 
200 
100 



HT 29 



Figure 4 | mRNA expression patterns of Katp channels in different cell Unes. (a) C3A; (b) A549; (c) A431; (d) HT29 cells. 



in the caveolae'", and more than one site in Kir6.2 (e.g., S372 or T224 
in Kir6.2)=' and SUR (e.g., T633 and S1465 in SUR2B)'''' have been 
impUcated in PKA phosphorylation. The positive regulation of K^tp 
by PKA was evidenced by the observed attenuation of Katp channel 
currents in the presence of PKA inhibitors'""^^. Of note, compounds 
containing 5-isoquinolinsulfonyl moieties such as H-89 also bind 
directly to the SUR subunit but with relatively low potency^\ At 
the concentration range tested it is not possible to distinguish 
between direct blocking of the Katp channel and indirect 
inhibition through abrogation of PKA activity. The inhibition of 
the pinacidU DMR by staurosporine may also be linked to PKA (Ki 
~ 7 nM), but the prevalence of other mechanisms cannot be ruled 
Nonetheless, these results suggest that the pinacidil DMR can 
be modulated by kinase activity. 

JAK3 and JAK2 are involved in Katp channel signaling. Janus 
kinases (JAKl, 2, 3) are protein tyrosine kinases involved in 
cytokine mediated cellular signaling and are crucial for a variety of 
cellular functions including cellular survival, proliferation, differen- 
tiation and apoptosis^'''". Given the attenuation of the pinacidil DMR 
by AG490, we examined the role of JAKs in Katp channel signaling 
using multiple assays. First, DMR inhibition assay showed that 
AG490 dose-dependently suppressed the pinacidil DMR, yielding a 
loglCso of -4.77 ± 0.10 (n = 3) (Fig. 7a). Second, siRNA knockdown 
studies showed that the treatment with two siRNAs for JAKl had 
little impact on the pinacidil DMR, but two siRNAs for JAK2 and two 
for JAK3 all markedly suppressed the pinacidil DMR (Fig. 7b and c). 
Third, western blotting showed that C3A cells primarily express 



JAK3, to less extend JAK2, but not JAKl. Moreover, the treatment 
of cells with 100 \xM pinacidil markedly increased the level of JAK2 
(Fig. 7d) and JAK3 (Fig. 7e). Importantly, the Kir6.2 siRNA 
treatment impaired the pinacidil induced increase of JAK3 protein 
level (Fig. 7f and g), suggesting that the pinacidil induced increase in 
JAK3 protein is due to the activation of Kir6.2 Katp- We also 
attempted, but failed, to detect p-JAK levels in the absence and 
presence of pinacidil. Several possible mechanisms are accounted 
for this; one may be due to the poor activity of anti-p-JAK 
antibodies used; another is that the pinacidil activated Katp did 
not cause phosphorylation of JAKs, but only increased the 
expression level of JAKs. Of note, we used moderately high dose 
(40 |J,M) of pinacidil for DMR profiling of the effect of kinase 
inhibitor and RNAi knockdown in order to improve sensitivity, 
but a higher dose (100 |iM) for JAK expression analysis to achieve 
maximal effect. Nonetheless, these results suggest the important role 
of JAK3, to less extent JAK2, in the Katp channel signaling in C3A 
cells. 

Katp channel signaling is Unked to ROCK and actin remodeling. 

Rho kinases (ROCKl and ROCK2) play important roles in the small 
GTPase RhoA initiated signaling pathways. Rho kinases are known 
to be involved in a variety of cellular functions including cytoskeleton 
organization, cell proliferation and apoptosis^**. Given the ability of 
Y27632 to attenuate the pinacidil response, we next examined the 
role of ROCK in the Katp signaling in C3A cells. DMR inhibition 
assay showed that Y27632 dose-dependently inhibited the pinacidil 
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Figure 5 | The RNAi knockdown effect on the DMR of 40 (iM pinacidil in C3A cells. (a,d) The DMR of pinacidil in untreated, mock transfection treated 
(mock) and scrambled RNAi (sc. RNAi) treated C3A cells, in comparison with the negative control, (b) The DMR of pinacidil in the mock transfection 
and three SURl RNAi treated cells, (c) The DMR of pinacidil in the mock transfection and three SUR2 RNAi treated cells, (e) The DMR of pinacidil 
in the mock transfection and three Kir6.1 RNAi treated cells, (f) The DMR of pinacidil in the mock transfection and three Kir6.2 RNAi treated cells. 
(g,h) statistical analysis of the effects of RNAi knockdown on the pinacidil DMR. *, p < 0.05; ***, p < 0.001. (i,j) The DMR of 32 nM pinacidil in C3A 
cells as a function of Katp blockers. The cells were first treated with each blocker at different doses for 1 hr, followed by stimulation with 32 |j.M pinacidil. 
(k) The impact of co-existence of 40 \M tolazamide on the potency of pinacidil. For (g-k), the pinacidil DMR amplitudes at 50 min post stimulation 
were plotted. Data represents mean ± s.d. (n = 6 for a-h; n = 3 for i-k). 



DMR, yielding a loglCso of -5.47 ± 0.07 (n = 3) (Fig. 8a). RNAi 
knockdown of ROCKl or ROCK2 was found to markedly attenuate 
the pinacidil DMR (Fig. 8b and c). Western blotting confirmed that 
the efficiency of ROCK RNAi knockdown was about 40-50% 
(Fig. 8d-f), consistent with the level of suppression of the pinacidil 
DMR by these siRNAs (comparing Fig. 8f with Fig. 8c) . Rho A activity 
results showed that the treatment of C3A cells with pinacidil had 



little impact on the RhoA activity (Fig. 8g), suggesting that the 
pinacidil activated Katp channels had no effect on the activity of 
RhoA, the upstream effector of ROCK. 

Given the important roles of ROCK signaling in actin remodeling, 
we finally examined the impact of microfilament modulators on the 
pinacidil response in C3A cells. Results showed that both cytochala- 
sin B and latrunculin A dose-dependently blocked the pinacidil 
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DMR, leading to a loglCso of - 5.87 ± 0.03 and - 5.8 1 ± 0.04 (n = 3) 
(Fig. 8h). Cytochalasin B acts by capping growing actin microfila- 
ments, while latrunculin A acts by sequestering actin monomers (Ki 
~ 0.2-0.4 |xM). In contrast, the microtubule disrupting agent noco- 
dazole up to 32 [iM had little effect on the pinacidil DMR (Fig. 8h). 
These results suggest that the K^tp channel initiated signaling in liver 
cells is linked to ROCK activity, leading to actin remodeling. 

Discussion 

We here present a strategy centered on label-free cell phenotypic 
profiling to discover the composition and signaling of an endogenous 
Katp channel in C3A cells. Label-free cell phenotypic profiling of a 
compound library consisting of seventy-two ion channel ligands 



across the five different cell lines led to identification of twenty- 
two ligands that displayed agonistic activity in at least one cell lines. 
Several classes of active ligands can be assigned based on their known 
mechanism of action (e.g., Ca^* mobilizing agents such as A23187, 
thapsigargin and cyclopiazonic acid), or using orthogonal assays 
(e.g., GPR35 agonists flufenamic acid, niflumic acid and NPPB). 

Combining DMR RNAi knockdown, ligand pharmacology 
profiling, and RT-PCR has allowed us to ascertain the presence, 
composition and signaling of a functional SUR2/Kir6.2 K^tp ion 
channel in C3A cells. Katp openers are a structurally diverse 
group of compounds, and also diverse in the modes of action. 
Diazoxide is a benzothiadiazine dioxide, pinacidil a cyanoguanidine, 
PCO-400 a benzopran, minoxidil a pyrimidine N-oxide sulfate, and 
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ZM226600 an anilide tertiary carbinol. The Kir6.2/SUR1 channel can 
be activated by diazoxide but is relatively insensitive to pinacidil, 
while the Kir6.2/SUR2A channel is only weakly responsive to diaz- 
oxide but can be activated by pinacidil^'. On the other hand, the 
Kir6.2/SUR2B channel can be activated by diazoxide™; ZM226600 
can activate a channel consisting of SUR2B co-associated with Kir6. 1 
or Kir6.2'''; minoxidil can specifically activate Kir6.1/SUR2 chan- 
nels". Our DMR profiling showed that among all KCOs tested in 
C3A cells pinacidil triggered a robust DMR, while PCO-400 and 
diazoxide are less active, ZM226600 and minoxidil were inactive. 
Considering gene expression, RNAi knockdown and K^tp inhibitor 
profiling data, our DMR results are better explained by the presence 
of a functional Kir6.2/SUR2A channel in C3A cells. Of note, due to 
the low expression of K^-pp channels, in particular Kir6.2 subunits at 
both mRNA and protein levels in C3A cells, we have difficulty to 
determine the exact location of these channels whose activation con- 
tributes to the DMR of pinacidil. We also did not determine the 
location and compositions of Katp channels in other pinacidil 
responsive cell lines. Nonetheless, these results demonstrate the 



use of label-free cell phenotypic profiling for broad pharmacological 
profiling of Katp channels. 

The molecular events underlying the modulation of Katp activity 
by protein kinases is poorly understood. Both JAK2 and ROCK 
inhibitors have been reported to be physiologically linked to Katp 
channels in cardioprotection. The JAK/STAT pathway mediates 
opioid induced cardioprotection via glycogen synthase kinase-3P; 
AG490 abrogates this cardioprotective effect*''. The proposed mech- 
anism for reduced myocardial no-reflow during ischemic precondi- 
tioning, a strategy also used for reperfusion therapy following acute 
myocardial infarction, involves the opening of Katp channels via 
inhibition of ROCKs*''. ROCK inhibitors decrease the area of no- 
reflow; Katp inhibitors abolish the reduced no-reflow. Although 
the exact cellular pathways undetermined, we here showed that both 
JAK3 (to less extent JAK2) and ROCK are involved in the Katp 
channel signaling in C3A cells, a liver cell line. 

Information regarding the role and pharmacology of Katp chan- 
nels in hepatocytes is limited. While the physiological outcomes of 
Katp channel activity in hepatocytes and cardiomyocytes are 
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expected to be different, it is likely that their modulation by kinases 
occur through similar molecular pathways. The ability to capture 
complex physiological phenomena surrounding ion channel activity 
in experimentally convenient cellular screening assays, as described 
herein, is powerful and unprecedented. Label-free DMR assays offer 
great flexibility in assay formats^^, rich information content with 
real-time kinetics, high throughput, and more importantly, holistic 
view of receptor signaling and drug pharmacology with wide cov- 
erage of pathway/targets^^, including ion channels as demonstrated 
here. DMR assays not only permit pharmacological profiling of dis- 
tinct types of ligands (openers and inhibitors) for the channel pro- 
teins, but also allow for investigating the signaling mediated through 
ion channels, later of which is difficult to be assessed using electro- 
physiology. However, DMR assays also come with several limita- 
tions. First, the label-free assays cannot detect directly the ion flux 
related to the open and close of ion channels. Second, the pharmaco- 
logical profiles of ligands for a specific ion channel may not correlate 
with those obtained using electrophysiology, given that the DMR is a 
whole cell phenotypic response linking to the functional responses of 
cells upon the activation of ion channels. Third, without validation 
using conventional approaches, label-free is difficult to differentiate 
the upstream (e.g., PKA phosphorylation of the Katp channels), and 
downstream (e.g., the Katp activation increased JAK expression) 
signaling events of the channels. Thus, label-free biosensor enabled 
cell phenotypic profiling for investigating ion channel pharmacology 
is an attractive complimentary technology to electrophysiology and 
other conventional techniques. 

Methods 

Materials. BioMol kinase inhibitor library and latrunculin A was obtained from Enzo 
Life Sciences, Inc. (Farmingdale, NY, USA). CID2745687, diazoxide, flufenamic acid, 
minoxidil, and pinacidil were obtained from Tocris Chemical Co. {St. Louis, MO, 
USA). PCO-400 was obtained from Santa Cruz Biotech Inc. (Santa Cruz, CA, USA). 
AG490, cytochalasin B, glibenclamide, glipizide, tolazamide, tolbutamide, U-27883A, 
and Y27632 were purchased from Sigma-Aldrich (St. Louis, MO, USA). CeU culture 
reagents were purchased from Life Technologies (Grand Island, NY, USA). The 
kinase inhibitor library, each at 10 mM in dimethyl sulfoxide (DMSO), were freshly 
diluted using the assay buffer (HBSS; 1 X Hanks' balanced salt buffer, 10 mM Hepes- 
KOH, pH 7.1) before use. All other compounds were dissolved in DMSO at 100 mM 
and stored at — SO^'C. Epic® 384-well biosensor cell culture compatible microplates 
were obtained from Corning Incorporated (Corning, NY, USA). 

Cell culture. AU five cell lines were obtained from American Type Cell Culture 
(ATCC) (Manassas, VA, USA). They are human epidermoid carcinoma A431, 
adenocarcinomic human alveolar basal epithelial ceU A549, human colorectal 
adenocarcinoma HT-29, human hepatocellular carcinoma HepG2 and its clonal 
derivative C3A (HepG2C3A). C3A is a clonal derivative of HepG2 that was selected 
for strong contact inhibition of growth, high albumin production, high production of 
alpha fetoprotein and ability to grow in glucose deficient medium. All cells were 
subcultured 1-2 times per week according to ATCC's instruction. Cell passage less 
than 15 was used for aU experiments. All cells were passaged at 37'^C with 5% CO2 
using the following complete medium: McCoy's 5A medium for HT-29, Eagle's 
Minimum Essential Medium (MEM) for HepG2 and HepG2C3A, F-12K Medium for 
A549, and Dulbecco's Modified Eagle's Medium (DMEM) for A43L All media were 
supplemented with 10% fatal bovine serum, 4.5 g/L glucose, 2 mM glutamine, 
100 |ig/ml penicillin and streptomycin. AIL cells were passed with trypsin/ethylene- 
diaminetetraacetic acid when approaching 90% confluence to provide new 
maintenance culture on T-75 flasks and experimental culture on the biosensor 
microplates. 

For DMR assays, the ceU seeding density was optimized for each cell line. The 
optimal seeding density was found to be 32 K,25 K,25 K,21 K,and21 Kper well for 
HT-29, A431, A549, HepG2 and C3A cells, respectively. After overnight culture 
(—20 hrs), aU cells were directly washed using the HBSS buffer before assay, except 
for A43 1 which was subject to starvation for an extra day using the serum-free DMEM 
medium. The cell confluency was —95% at the time of assaying for all cells. 

Reverse transcriptase PGR. Total RNA was extracted from one T-75 flask with a 
confluent monolayer of cells (15-30 million cells) for each cell line using an RNeasy 
mini kit (Qiagen, Cat#74104, Valencia, CA). On-column DNase digestion was 
performed using RNase-free DNase set (Qiagen, Cat#79254) to eliminate genomic 
DNA contamination. The concentration and quality of total RNA were determined 
using a Nanodrop 8000 (Thermo Scientific). The primer sequences used were based 
on previous reports^^'^^ and custom synthesized through Sigma-Aldrich. Reverse 
transcriptase PGR was performed using One-Step RT-PGR kit from Qiagen Inc. 
(Cat# 210212). The PGR conditions were as follows: 50 G for 30 min, 95 "C for 



15 min, followed by 60 cycles of 1 min at 94''G, 1.5 min at 57 G and 2 min at 72°G, 
with a final extension of 10 min at 72 'C. For aU gel electrophoresis analysis equal 
amount of DNA was loaded. The DNA of actin obtained was used the control. The 
comparative level was estimated across different genes in a specific cell line, not across 
different cell lines. 

RNAi knockdown. All siRNAs including universal scrambled RNAi were selected 
from the pre-designed siRNA database of Sigma-Aldrich. For each target the top two 
or three siRNA having predicted high knockdown efficiency was selected. The siRNA 
IDs are SASI_Hs01_00092347, SASI_Hs01_00092348 and SASI_Hs01_00092349 for 
SURl (Refseq ID: NM_000352); SASI_Hs01_00061301, SASI_Hs01_00061302 and 
SASI_Hs01_ 00061303 for SUR2 (Refseq ID: NM_005691); SASI_Hs01_ 00113357, 
SASI_Hs01_001 13358 and SASI_Hs01_001 13359 for Kir6.1 (Refseq ID: NM004982); 
SASI_Hs01_00220256, SASI_ Hs01_00220257 and SASI_Hs01_00220258 for Kir6.2 
(Refseq ID: NM000525); SASI_Hs01_00065573, SASI_Hs01_00065571, and 
SASI_Hs01_00065570 for ROGKl (RefseqID: NM_005406); SASI_Hs01_00204253, 
SASI_Hs01_00204252, and SASI_Hs01_00204251 for ROCK2 (RefseqID: 
NM_004850); SASI_Hs01_ 00174614 and SASI_Hs01_ 00174613 for JAKl (ReqID: 
NM_002227). SASI_Hs02_ 00338675 and SASI_Hs01_00041547 for JAK2 (ReqID: 
NM_004972); SASI_Hs01_001 18128 and SASI_Hs02_ 00302103 for JAK3 (ReqID: 
NM_000215). 

siRNA transfection was performed using the N-TER nanoparticle siRNA 
Transfection system (Gat# N2913) from Sigma. Specifically, 5000 cells were first 
plated into each well of an Epic® 384well microplate, and cultured for 20 hrs using the 
complete medium. The cells were then transfected with 50 nM siRNA and incubated 
in siRNA containing medium for 24 hrs before replaced with fresh ceU culture 
medium. Label-free DMR assays were performed 48 hrs after transfection. The cells 
treated with the transfection vehicle were used as the mock control. 

Immunoprecipitation and Western blotting. G3A cells were plated in 6-weIl tissue 
culture treated plate with 3 X 10^ cells per well. Cells were transfected with respective 
top-rank siRNA at 50 nM final concentration 20 hrs after plating. 24 hrs after 
transfection the transfection reagent containing media were removed and replaced 
with the respective complete media. 48 hrs after transfection cells were lysed using 1% 
NP40 lysis buffer (150 mM NaCl, 25 mM Tris-HGl, pH 7.6, 1% NP-40) with the 
protease inhibitor cocktail (Roche Applied Science, Indianapolis, IN, USA). The cell 
lysate was then centrifuged at 14,000 rpm for 20 min. The supernatant was then 
transferred to a new tube before western blotting. 

For west blotting, 120 |il cell lysate of each sample were mixed with 40 |,d 4X SDS 
sample buffer (Bio-Rad Life Sciences, Hercules, CA, USA), and then boiled at lOO^C 
for 5 min. Proteins were separated on 10% SDS gel, and 15 |j.l of each sample was 
loaded to the gel. For ROGKl and 2, membrane was blotted with either rabbit anti- 
ROCKl (Santa Cruz Biotech, catalog # sc-5560), or rabbit anti-ROCK2 (sc-5561), or 
goat anti- Actin (sc-1616) (1 : 500 dilution) for 1 hr, and then with 2nd HRP con- 
jugated goat anti-rabbit or horse anti-goat antibody (1 :2000 dilution) for 15 min. 
Similarly, for JAKs, cell lysate was immunoprecipitated with rabbit anti-JAKl (sc- 
7228), goat anti-p-JAKl (Tyr 1022/Tyr 1023) (sc-16773), rabbit anti-JAK2 (sc-294), 
goat anti-p-JAK2 (Tyr 1007/Tyr 1008) (sc-21870), rabbit anti-JAK3 (sc-513) or goat 
anti-p-JAK3 (Tyr 980) (sc- 16567) (1 : 250) at 4"G and then with 2nd HRP conjugated 
goat anti-rabbit or horse anti-goat antibody (1 : 2000 dilution). The images were 
developed with ECL-Plus (GE Healthcare). The original west blot images were pre- 
sented in Supplementary Figures SI to S5. 

DMR assays. DMR assays in microplate were performed using high throughput 
screening compatible Epic® system (Corning)^^. This system consists of a 
temperature-control unit (26'^G), an optical detection unit, and an on-board liquid 
handling unit operated by robotics. This system scans entire biosensor microplate 
every 6 sec, and two scans are averaged to reduce the signal noise, thus leading to a 
kinetic response with a temporal resolution of — 15 sec. The cultured cells were first 
washed with the assay buffer HBSS using a plate washer (Bio-Tek Microplate Washers 
ELx405"^, Bio-Tek, Winooski, VT), and incubated inside the Epic® system for 
— 1 hr. After a steady baseline was established, the ligand solutions were introduced 
using the on-board liquid handling device, and the cell responses were then recorded 
over time. For inhibitor profiling, the cells were pretreated with an inhibitor for 1 hr 
before pinacidil stimulation. For RNAi, the cells were transfected with siRNA for 
48 hrs before pinacidil stimulation. All real-time DMR signals were reported as a 
2 min baseline, followed by ligand stimulation. All DMR signals were background 
corrected. At least two independent sets of experiments, each with at least duphcate, 
were performed. 

Tango p-arrestin translocation gene reporter assay. Tango assays were performed 
in engineered Tango™ U20S-GPR35-bla cell line (Life Technologies). This cell line 
allows an endpoint measure of the activity of agonists specific to the GPR35 
activation -induced P-arrestin translocation^^-™. The cells were passed using McCoy's 
5A medium supplemented with 10% dialyzed fetal bovine serum, 0.1 jiM NEAA, 
25 )j.M Hepes (pH 7.3), 1 mM sodium pyruvate, 100 U/ml penicillin, 100 jJ-g/ml 
streptomycin, 200 (ig/mlzeocin, 50 jig/mlhygromycin, and 100 |jg/ml geneticin in a 
humidified 37' C/5% CO2 incubator. For Tango assays, 10000 \}10^-Gm35-h\a cells 
per well were seeded in 384-well, black-wall, clear bottom assay plates with low 
fluorescence background (Corning). After overnight culture, the cells were stimulated 
with ligands for 5 hrs at 37^C under 5% GO2, and then loaded with the cell permeable 
LiveBLAzer™ FRET (fluorescence resonance energy transfer) B/G substrate. After 
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2 hr incubation the coumarin to fluorescein ratio was measured using Victor 4 plate 
reader (PerkinElmer, Waltham, MA, USA). Results obtained were normalized to the 
maximal response of zaprinast {set to be 100%). 

RhoA activation assay. RhoA activity was determined from protein isolated from 
C3A cells without and with treatment with 100 |j.M pinacidil using the luminescence 
based G-LISATM RhoA activation assay kit (Cytoskeleton, Denver, CO, USA) 
according to the manufacturer's instructions. Protein was isolated using the provided 
cell lysis buffer, and cells were processed rapidly on ice and snap-frozen until the time 
of assaying. Lysates were clarified by centrifugation at 10,000 rpm at 4'"C for 2 min. 
Protein concentration was determined according to the manufacturer's protocol, and 
cell extracts were equalized to contain protein concentrations of 2 mg/ml for the 
assay. Luminescence was detected as suggested by the manufacturer. 

Data analysis. All dose- dependent responses were analyzed by using GraphPad 
Prism 5.0 (GraphPad Software Inc., San Diego, CA, USA). The EC50 and IC50 values 
were obtained by fitting the DMR dose response curves with nonlinear regression. For 
cluster analysis, we profiled the ion channel ligand library across the five cell lines. For 
effective clustering, we extracted the real responses at the six distinct time points {3, 5, 
9, 15, 30, 45 min) for each DMR to form a numerical description of the label-free cell 
phenotypic characteristics^^ All the six time points refer to the stimulation 
duration after renormalizing the responses starting from the time when the 
compound was added. The real responses at these discrete time points were color 
coded to illustrate relative differences in the direction and strength of a DMR signal. 
The red color indicates a positive value, while the green color refers to a negative 
value, the black a value at or near zero. Differences in color intensity illustrate 
differences in signal strength. In the ligand heatmap matrix (Fig. 1) each column 
represents one DMR value at a specific time in a specific cell line, and each row 
represents one ligand. The Ward hierarchical clustering algorithm and Euclidean 
distance metrics (http://www.eisenlab.org/eisen/) were used for similarity analysis, 
and every row and column carries equal weight. DMSO in the assay buffer at a 
concentration that equals to those for all ligands was also included as a negative 
control. Each DMR represents an average of four replicates. To assist with direct 
visualization of DMR characteristics of each ligand in a cell line, we did not carry out 
similarity analysis among distinct columns. 

Statistical analysis. DMR data were analyzed by using GraphPad Prism 5.0 
(GraphPad Software Inc., San Diego, CA, USA). The EC5Q values were obtained by 
fitting the dose DMR response curves with nonlinear regression. The RNAi effect was 
analyzed using one-way ANOVA Tukey analysis with the Prism 5.0. 
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